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1.  INTRODUCTION 


The  spontaneous  penetration  of  liquids  into  porous  materials  under  the  influence  of  capillary  forces 
is  of  importance  in  many  physical  and  physiological  processes.  Understanding  the  basic  mechanisms  and 
relevant  parameters  associated  with  impregnation  is  required  in  fields  as  diverse  as  materials  science  and 
geology.  Because  of  their  widespread  use  in  many  advanced  composite  materials,  we  examine  the 
impregnation  of  tows  of  carbon  fibers.  Our  particular  interest  relates  to  the  manufacturing  of  continuous 
fiber  composite  materials,  wherein  an  essential  unit  process  consists  of  impregnating  bundles  of  dry  fibers 
(a  fiber  tow)  with  a  viscous  resin.  A  typical  carbon  fiber  tow  is  composed  of  12,000  individual  fibers  (a 
12K  tow)  having  perimeters  from  15  to  30  pm.  The  circumferences  of  the  fibers  we  investigate  are  nearly 
circular.  The  tows  constitute  a  simple  model  porous  system  of  axially  aligned  rods  and  we  consider 
impregnation  driven  solely  by  capillary  forces. 

In  a  regular  array  of  parallel  fibers,  two  types  of  flow  can  be  envisaged:  axial  flow  along  the  fiber 
axis  and  transverse  (normal)  flow  normal  to  the  axes  of  the  fibers.  In  a  perfectly  regular  lattice  structure 
of  fibers  (e.g.,  hexagonaUy  packed),  axial  penetration  occurs  if  the  fibers  touch  with  their  axes  normal  to 
the  liquid  surface.  Normal  impregnation  takes  place  if  the  fibers  touch  with  their  axes  parallel  to  the 
liquid  surface.  In  many  cases  of  practical  importance  both  modes  of  impregnation  co-exist  due  to  (1)  the 
irregularities  in  the  lattice  structure  and  (2)  the  application  of  liquid  in  an  arbitrary  direction.  However, 
in  processes  associated  with  the  manufacture  of  continuous  fiber  composites,  the  normal  impregnation 
mode  dominates.  The  adhesive  is  applied  to  the  surface  of  a  continuous  tow  and  impregnation  occurs  by 
a  combination  of  capillary  forces  and  applied  pressure. 

In  a  previous  publication  (Bayramli  and  Powell  1990),  we  studied  normal  impregnation  theoretically 
and  predicted  rates  10  to  100  times  slower  compared  with  axial  impregnation  as  described  by  the  Lucas- 
Washbum  approximation  (Lucas  1918;  Washburn  1921).  An  experimental  study  of  normal  impregnation 
has  provided  qualitative  confirmation  of  these  theoretical  predictions  (Bayramli  and  Powell,  to  be 
published). 

In  this  study,  we  investigate  axial  impregnation  experimentally  and  compare  the  results  with  the 
theoretically  predicted  rates  of  wicking  in  a  hexagonally  packed  regular  array  of  fibers.  For  all  data 
reported,  the  capiUary  number  (ratio  of  the  surface  to  gravitational  forces)  is  sufficiently  large  (>  10^)  so 
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that  surface  tension  forces  dominate.  A  novel  analysis  of  the  data  permits  the  calculation  of  the  effective 
pore  size  and  the  pore  size  distribution  in  the  fiber  tows. 

2.  THEORETICAL  BACKGROUND 

We  investigate  flow  in  a  porous  medium  which  has  a  constant  cross-sectional  geometry  normal  to  the 
direction  of  flow.  Figure  1  shows  such  an  array,  with  the  flow  occurring  normal  to  the  plarre  of  the  figure 
(i.e.,  along  the  fibers’  axes).  It  is  assumed  that  the  fluid  is  linearly  viscous  (with  viscosity  4),  that  inertial 
and  gravitational  effects  can  be  neglected,  and  that  the  flow  is  driven  solely  by  capillary  forces.  If  at 
r  =  0  impregnation  commences,  then  the  distance  of  impregnation  of  the  liquid  into  the  medium,  h,  is 
given  by  Lucas  (1918)  and  Washburn  (1921), 

h  =  (1) 

A=A'(y  cos  e/4/"  (2) 

where  A'  is  a  constant  which  basically  depends  upon  the  geometry  of  the  medium,  y  is  the  surface  tension, 
and  9  is  the  parent  contact  angle  of  the  contact  line  on  the  solid  waU. 

For  simple  models  of  porosity,  it  is  possible  to  determine  the  functional  form  of  A\  For  example,  if 
the  interstices  of  the  medium  consist  of  uniform  hollow  cylinders  with  radius  r,  then 

/l'  =  (r/2)'"  .  (3) 

If  the  gemnetry  is  more  complicated  (e.g.,  a  regular  hexagonal  packing  of  axially  oriented  solid  cylinders, 
where  the  flow  occurs  in  the  pores  between  the  cylinders),  the  calculations  are  cumbersome.  Not  only 
is  tire  pressure  jump  across  the  interface  hard  to  calculate  because  of  its  complicated  shape,  but  also  the 
solution  of  the  creeping  flow  equations  is  difflcult  due  to  the  geometrically  complex  boundary  conditions. 
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Note:  The  r,  and  r,  denote  the  Radii  of  two  imaginary  cainllaries  and  2d  is  the  minimum  fiber-to-fiber  spacing.  Parallel 
plates  extend  along  the  flber’s  axis  and  the  body  angle  $  is  increased  in  such  a  way  as  to  compensate  one-third 
of  the  volume  of  the  shaded  area  (Bayramli  and  Powell  1990).  The  other  two-thirds  is  compensated  by  the  neighboring 
two  pares. 


Figure  1.  Regular  Hexagonal  Packing  of  Fiber. 


It  is  possible  to  calculate  the  constant  ^4'  as  a  function  of  porosity  for  various  pore  geometry 
approximations  of  the  regular  hexagonal  packing  of  fibers  shown  in  Figure  1 .  The  porosity  P  is  calculated 
for  the  unit  triangular  cell  DEF  in  Figure  1  as 


(2v/3)  [1  +  (d/R)]" 


(4a) 


where  R  is  the  fiber  radius  and  2d  is  the  minimum  fiber-to-fiber  separation. 

If  the  impregnation  dynamics  are  approximated  by  assuming  that  wicking  occurs  in  a  series  of  hollow 
cylinders  between  the  fibers,  impregnation  rates  A,  and  Ai  corresponding  to  and  rj  (Hgure  1)  can  be 
calculated.  For  diese  simple  geometries,  the  imaginary  cylindrical  capillaries  with  radii  r,  and  rj  can  be 
related  to  /*  as  follows; 


[rt/?V2^(l  -  P)Y^ 
cos;c/6 


(4b) 


3 


tan— 
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JC/J* 


2^3  (1  -  P) 


\in 


(4c) 


The  rate  A  is  then  calculated  as  a  function  of  porosity  using  Equations  4a,  4b,  and  4c  and  substituting  into 
Equations  1  and  3,  noting  that  both  and  rj  are  uniquely  determined  by  d. 

It  is  most  useful  to  express  our  experimental  results  as  well  as  the  model  predictions  in  dimensionless 
terms.  Expressed  nondimensionaUy,  Equation  1  is 


h 


^  At 


1/2 


(5) 


where  the  bars  denote  nondimensional  quantities  which  are  defined  as 


h 


-  ^  YCOS0 , 

"'~W~ 


The  parameter  A  can  be  calculated  using  a  model  which  divides  the  fiber  perimeter  into  parallel  plates 
extending  along  the  fiber  axis,  as  shown  in  Figure  1.  By  this  method,  we  can  make  use  of  the  actual 
capillary  pressure  for  hexagonal  packing  derived  by  Princen  (1969),  rather  than  using  the  cylindrical  pore 
approximation.  The  viscous  resistance  is  calculated  by  dividing  the  fibers’  surfaces  into  small  parallel 
plates  extending  along  the  axes  of  the  fibers.  The  excess  shaded  triangular  area  in  Figure  1  can  be 
incoqporated  into  the  parallel  plate  approximation  by  inc.xasing  the  body  angle  ({>  slightly  so  that  all  the 
pore  area  is  takoi  into  account  For  pressure  driven  flow  along  the  parallel  plates,  the  flow  rate  depends 
on  the  cube  of  the  separation  of  the  plates,  2w.  Therefore  a  flow  average  separation,  2w^  can  be  defined 
as 


2w 


1 

n 


m 


|U3 


E  (2h'.)' 


(6) 
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where  /I  is  the  number  of  parallel  channels  which  follow  the  contour  of  die  fiber  perimeter  closely  (an 
arbitrarily  large  number  chosen  to  be  1,(XX)).  For  the  triangular  unit  cell  in  Figure  2,  DEF,  Princen  (1969) 
defined  the  upward  surface  tension  force  for  porosities  above  0.20  for  6  =  0°  C  and  above  0.153  for 
6  =  60°.  From  his  results,  the  capillary  pressure  (^)  for  the  hexagonal  geometry  can  be  easily  obtained 
by  dividing  it  by  the  pore  area  over  which  it  acts: 

(7) 

R[(l  +  -  jt/2] 


where  d  =  d/R .  Substituting  Ap  from  Equation  7  into  the  Lucas-Washbum  type  creeping  flow  equation 
for  the  parallel  plate  approximation  of  the  hexagonal  fiber  packing,  we  And 

1/2 

(8) 


2w. 


1C  cos  6 


6(1  +df>l3  -  Jt/2 


Here  and  is  the  impregnation  constant  for  the  hexagonally  packed  regular  fiber  structure 

appearing  in  an  appropriately  dimensionless  form  of  Equation  S. 

The  theoretical  impregnation  rates  for  fiber  bundles  having  various  porosities  as  a  function  of 
i4, ,  /4j,  are  compared  with  the  experimental  results  in  Section  4. 

3.  EXPERIMENTAL 

3.1  Materials.  The  mass  of  liquid  impregnated  into  axially  oriented  carbon  fiber  tows  is  measured 
by  suspending  the  tow  from  an  electrobalance  with  the  fibers’  axes  aligned  with  gravity.  The 
electrobalance  is  a  Calm  RlOO  which  can  measure  to  1  pg  and  is  stable  at  10  pg  up  to  one  day.  A  beaker 
containing  the  liquid  is  placed  on  a  movable  stage.  The  beaker  is  raised  until  the  surface  of  the  liquid 
touches  the  tow  and  forms  a  meniscus  around  it  (Figure  2a).  All  experiments  reported  here  were 
conducted  at  room  temperature.  The  carbon  fiber  samples  for  the  impregnation  studies  are  prepared  by 
placing  a  0.4-m'long  tow  under  tension  and  applying  adhesive  tape  tightly  to  both  sides  of  the  tow  in  a 
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(a) 


(b) 


Note:  (a)  Tow  impregnation  and  single-fiber  experiments,  and  (b)  the  preparation  of  tow  samples  under  tension.  W  is 
0.239  kg. 

Figure  2.  Schematic  Diagram  of  the  Experimental  Setup. 

way  that  encapsulates  it  (Figure  2b).  The  strip  is  then  cut  into  pieces  20-40  mm  long,  and  each  of  these 
is  trimmed  on  the  sides  to  be  used  in  an  impregnation  experiment  The  taping  process  ptxxluces  a 
controlled  volume  and  gives  rigidity  to  the  sample.  It  can  be  suspended  from  one  of  the  anns  of  the 
electrobalance  through  a  small  hole  drilled  in  the  sample.  The  liquid  enters  the  sample  only  through  the 
bottom  cross-sectional  surface  located  at  the  hydrostatic  zero  level.  A  blank  test  made  by  adhering  two 
tapes  face-to-face  and  trimming  down  to  sample  size  showed  no  impregnation.  Further,  the  interference 
of  the  adhesive  tape  with  the  surface  tension  measurement  can  be  neglected,  and  the  contact  angle  of  the 
silicone  oil  on  carbon  fiber  can  be  taken  to  be  zero.  For  the  epoxy  system,  a  certain  amount  of 
contamination  is  possible  through  dissolution  of  the  adhesive  on  the  tape.  However,  no  change  in  the 
measured  wetting  force  of  the  epoxy  system  was  observed  when  it  was  in  contact  with  the  adhesive  tape. 
Lastly,  the  ratio  of  the  carbon  fiber  surface  area  inside  the  sample  compared  with  the  surface  area  of  the 
enclosing  tape  is  about  25/1.  Therefore,  depending  on  the  contact  angle  of  the  liquid  on  die  tape,  it  can 
introduce  a  maximum  of  4%  error  in  the  measurements. 

The  carbon  fibers  are  availaUe  as  continuous  12K  tows.  The  particular  fibers  and  their  suppliers  are 
listed  in  TaUe  1.  All  results  reported  here  are  for  sized  fibers.  The  perimeters  of  single  carbon  fibers  are 
measured  by  the  fiber  balance  wetting  technique,  which  has  now  become  a  standard  way  to  characterize 
single  fibers  (Okagawa  and  Mason  1978).  For  this  purpose,  hexadecane  is  used  as  die  wetting  liquid, 
which  gives  a  neariy  zero  contact  angle  with  carbon  fibers.  The  experimental  setup  is  essentially  the  same 
as  that  used  for  the  impregnation  studies  and  is  shown  schematically  in  Figure  2.  Contact  angles  of  the 
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Table  1.  Properties  of  the  Fibers  and  Liquids 


CARBON  FIBER 

Perimeter* 

Density 

(pm) 

(kg  m-') 

IM7  (Hercules) 

15.7 

1,870 

21.6 

Magnamite  AS4  (Hercules) 

25.1 

1,770 

21.2 

1,780 

LIQUID 

Viscosity 

Density 

Surface  Tension 

(Pas) 

(kg  m-^) 

(mN  m’*) 

Silicone  oil*’  (10  cs) 

0.102 

940 

20.1 

Silicone  oD**  (100  cs) 

0.93 

968 

20.6 

Silicone  oil”  (500  cs) 

4.75 

972 

21.1 

Hexadecane 

773 

25.7 

Adhesive”^ 

13.4 

1,102 

35.3 

*  Average  value. 

^  Dow  Coming  200  fluids;  values  at  25°  C. 

°  100/40  (w/w)  Epoxy  (bisphenol-A-di-glycidyl  ether)/hardener  (aliphatic  polyethertriamine). 

**  Viscosi^,  density,  and  surface  tension  of  the  epoxy  system  are  about  constant  for  3  hotirs  after  mixing. 


epoxy  system  with  the  carbon  fibers  are  determined  from  the  perimeter  measurements  and  the  total  wetting 
force  on  single  fibers. 

The  densities  of  the  carbon  fibers  are  determined  using  a  pycnometer.  A  1-m-long  tow  is  cut  into 
2-mm-long  strips  and  completely  wetted  in  hexadecane.  From  the  weight  of  the  pycnometer  with 
hexadecane,  and  with  hexadecane  plus  carbon  fiber,  the  density  of  the  carbon  fibers  can  be  calculated. 
The  weight  per  unit  length  of  a  tow  is  measured  by  applying  the  same  tension  as  that  used  in  the 
preparation  samples  for  the  impregnation  studies  to  a  nominally  l-m-long  tow  and  measuring  the  exact 
length  and  weight  of  the  tow. 

The  epoxy  system  is  bisphenol-A-di-glycidyl  ether  (Der  332,  Dow  Chemicals)  and  an  aliphatic 
polyethertriamine  (Jeffamine  T403,  Jefferson)  at  l(X)/40  (w/w)  epoxy/hardener  ratio.  It  is  degassed  at 
0.35  atm  after  mixing,  and  a  profuse  amount  of  bubble  formation  is  observed  which  lasts  about  15  min. 
In  addition  to  the  epoxy,  various  silicone  fluids  were  used  (Table  1).  The  surface  tension  of  the  liquids 
used  in  the  experiments  is  determined  by  the  Wilhelmy  Plate  technique  and  the  viscosities  are  measured 
using  a  Visco-electric  Analyzer  (Sangamo-Schlumberger,  Bognor  Regis,  U.K.). 
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Cross-sectional  views  of  the  epoxy  impregnated  aiul  cured  tows  are  obtained  by  optical  reflectance 
and  scanning  electron  microscopies.  The  surfaces  of  the  samples  are  prepared  by  diamond  polishing. 

3.2  Analysis.  The  result  of  a  typical  experiment  is  given  in  Figure  3  where  the  nondimensional 
volume  V  •  V/R^  (V  being  the  instantaneous  volume  of  liquid  impregnated)  is  shown  as  a  function  of 
the  square  toot  of  the  dimensionless  time  t .  In  this  and  subsequent  figures,  V  is  calculated  by  dividing 
the  measured  mass  by  density.  The  volume  increases  lineaily  with  7  ^ .  Towards  the  end  of  the 
experiment,  there  is  deviation  from  linearity,  indicating  that  the  impregnation  front  in  the  tow  is  irregular. 
Capillaries  with  different  pore  dimensions  have  different  wicking  rates.  There  is  attrition  as  those  with 
the  large  diameters  (with  higher  rates  of  impregnation)  are  completely  filled  while  flow  continues  in  pores 
of  smaller  dimensions.  The  net  effect  is  to  decrease  the  overall  rate  of  impregnation.  We  note  that 
nonlinear  behavior  may  result  also  from  levelling  off  as  the  equilibrium  capillary  height  is  approached. 
In  these  experiments,  the  samples  are  sufficiently  short  so  that  such  effects  are  not  present 

This  finding  is  further  supported  by  the  studies  using  epoxy.  The  epoxy  system  used  in  these 
experiments  has  a  long  pot  life  with  its  room  temperature  viscosity  staying  practically  unchanged  for 
3  hours  after  mixing  (Morgan,  Kong,  and  Walkup  1984).  Before  it  hardens  completely,  the  maximum 
distance  the  curing  epoxy  system  can  axially  impregnate  into  a  carbon  fiber  tow  is  about  40  mm, 
depending  on  the  type  and  the  porosity  of  the  tow.  An  examination  of  the  hardened  cross  sections  of 
epoxy- impregnated  tows  near  the  leading  edge  of  the  meniscus  reveals  a  penetration  difference  of  as  much 
as  4  mm  in  various  pores  of  the  samples.  It  is,  therefore,  not  possible  to  relate  accurately  the  volume 
increase  in  the  linear  part  of  figure  3  to  a  unique  contact  line  displacement.  However,  an  estimate  of  the 
overall  wicking  rate  can  be  obtained  from  the  slope  of  the  linear  region  of  the  impregnation  curve. 

It  is  not  possible  to  detennine  the  initial  point  on  the  curve  in  Figure  3  by  subtracting  the  external 
meniscus  force  from  the  total  measured  force.  This  results  from  the  meniscus  force  being  large  with 
respect  to  the  weight  of  the  liquid  inside  the  tow  at  early  times,  and  from  the  inaccuracy  in  the 
determination  of  the  external  perimeter  of  the  tow  sample.  The  initial  point  shown  in  Figure  3 
(7  *  0,  V  >  0)  is  calculated  by  extrapolating  the  measured  volume,  after  the  formation  of  the  meniscus, 
back  to  r  =  0  assuming  a  linear  V  vs.  t  relationship.  If  the  linearity  assumption  were  not  valid,  the 
straight  line  would  not  pass  through  subsequent  data  points,  which  is  found  to  be  the  case  in  all 
experiments. 
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Note:  Using  Equation  9,  the  porosity  is  calculated  to  be  0.S88.  The  slopes  at  points  1  and  2  are  for  illustration  purposes  and 
show  Vf,  from  which  the  pore  size  distribution  is  calculated.  The  maximum  amount  of  liquid  volume  impregnated  is  V^. 

Figure  3.  Impregnation  Curve  of  a  40-cm-Long  IM7  Carbon  Tow  With  100  cs  Silicone  Oil. 

The  porosity  of  the  tow  can  be  calculated  from  the  maximum  volume  of  liquid  impregnated  (V^),  the 
length  of  the  tow  (L),  the  density  of  the  carbon  fibers  (pp),  and  the  mass  per  unit  length  of  the  tow  (Mj) 
using 


Pf 


(9) 


The  experimental  rate  of  impregnation  can  be  calculated  approximately  from  the  linear  portion  of 
Figure  3,  assuming  the  front  advances  smoothly  with  NL  =  VIV^.  A  more  accurate  experimental  rate  is 
calculated  by  making  use  of  the  nonlinear  portion  of  Figure  3.  If  the  lateral  motion  of  the  liquid  in  the 
tow  is  ignored,  the  fiber  bundle  can  be  visualized  as  a  collection  of  individual  capillaries  with  a 
distribution  of  pore  sizes.  The  volume  of  liquid  impregnated  to  a  maximum  distance  within  pores  that 
have  characteristic  dimensions  between  r  and  r  +  Ar  is  determined  by  extrapolating  the  tangent  of  the 
V  vs.  curve  back  7*0. 

To  prove  this,  let  N(y)  denote  the  number  of  pores  having  volume  V,  and  be  the  largest  pore 
volume.  As  time  progresses,  the  larger  pores  tend  to  fill  first  The  total  volume  of  liquid  in  a  sample, 
V^,  is  given  by 
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v,(t) 


V^it)  -  J  N^)dV  +  J  N(y)dV 


V'lCO 


(10) 


In  Equation  10,  die  secraid  tenn  accounts  for  the  volume  in  those  capillaries  which  have  completely  filled 
at  time  t,  while  the  volume  contribution  from  the  first  term  is  from  capillaries  in  which  lit^d  is  still 
flowing,  with  Vj(r)  being  the  largest  capillary  in  whidi  there  is  flow.  We  consider  our  volumes  as  being 
composed  of  a  series  of  capillaries,  for  which  by  Equations  1-3,  when  flow  occurs 

V  =  (11) 

dV  ^  .  (12) 


Substituting  Equation  12  into  Equation  10, 


r^(t) 

J  /V (7)7 dr 
.  0 


+  J  N(r)h^2nrd'r 

7,(0 


(13) 


where  7,  (r)  is  the  largest  pore  radius  being  filled  at  time  t ,  and  7^  is  the  maximum  pore  radius  in  a 
sample.  Differentiating  Equation  13  with  respect  to  r  and  ^plying  the  Leibniz  rule 


dVj. 


7,(0 

J  N(r)7^dr 


(14) 


Equation  14  is  the  tangent  to  the  V,.  vs.  t  curve.  Extrapolating  the  tangent  back  to  r  =  0  and  letting  V  (r) 
be  tlw  intercept,  a  comparison  of  the  integrated  form  of  Equation  13  and  Equation  14  gives 
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(15) 


This  is  simply  the  contribution  of  completely  filled  capillaries  to  the  overall  volume  at  time  t . 


To  calculate  a  pore  size  distribution,  we  consider  two  intercepts.  V.  _ ,  and  V.,  with  A  V.  =  V.  -  V._^. 
Letting  N.  -  Nir.)  be  the  number  of  pores  having  volumes  between  V . .  j  and  V.,  and  employing 
Equation  15.  we  have 


IV.-  = 


(16) 


where  r.  is  the  mean  radius  of  the  capillaries  having  volume  AV..  To  use  Equation  16.  we  begin  by 
calculating  N(r.)  at  the  time  when  all  the  capillaries  are  full.  At  regular  7  intervals,  slopes  arc 
taken  and  extrapolated  to  «  0.  The  first  non-zero  slope  belongs  to  an  average  of  the  smallest 
capillaries  when  all  the  other  larger  capillaries  arc  full;  therefore,  it  represents  a  rate  which  is  solely  due 
to  the  smallest  capillaries.  The  first  (V^  -  V,)  gives  the  volume  of  these  small  pores.  Subsequent  slopes 
and  the  volume  increments  belong  to  progressively  larger  pores  and.  hence,  independent  contributions 
(V,  -  V,  _ ;)  of  all  individual  mean  capillaries  to  the  maximum  volume.  V^,,  can  be  determined.  The 
number  of  capillaries  can  then  be  calculated  from  Equation  16  employing  the  rate  relationships  given  in 
Equations  1-3  for  the  case  of  cylindrical  geometry.  For  an  accurate  extrapolation  of  the  tangents  in  the 
nonlinear  part  of  the  curve  back  to  the  V  axis,  a  smoothing  of  the  data  is  performed  using  a  polynomial 
approximation  of  the  curve.  It  is  found  that  a  fourth-order  polynmninal  fits  the  nonlinear  data  points  and, 
therefore,  an  accurate  extrapolation  of  the  data  points  to  the  volume  axis  can  be  made.  We  have  divided 
the  nonlinear  portion  of  the  curve  from  B  to  C  into  SO  equal  7  increments  and  obtained  SO  extrapolated 
values  Of  Vf  such  that  the  total  volume  impregnated  at  long  times.  Vj.(«>),  is  given  by 
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(17) 


Vr(^)  -  E  AV,  . 


■  >  1 


The  entire  curve  pan  in  Figure  3  from  B  to  C  is  used  to  calculate  an  average  experimental 
impregnation  constant  A  as  follows: 


so 


i  «  1 


(18) 


It  should  be  noted  that  no  assumption  as  to  the  pore  geometry  is  required  to  find  A  experimentally,  but 
to  compare  with  theory,  a  model  pore  geometry  must  be  assumed. 

4.  RESULTS  AND  DISCUSSION 

Among  the  carbon  fibers  used  the  IM7  has  the  most  tegular  and  round  cross  section  (Figure  4).  We 
have  used  silicone  liquids,  which  have  contact  angles  of  approximately  zero  on  carbon  fibers,  to  determine 
the  kinetics  of  axial  impregnation  into  IM7  fiber  bundles.  Representative  results  are  summarized  in 
Table  2.  The  porosities  in  experiments  1-9  varied  from  49.65  to  59.60.  The  experimental  values  of  A 
calculated  from  the  linear  and  the  nonlinear  portions  of  the  impregnation  curve  differed  slightly. 
Therefore,  a  good  approximation  can  be  achieved  by  using  h  vs.  t  of  the  linear  part  in  the  kinetic 
relationships.  The  data  in  Table  2  demonstrate  that  A  is  tx)t  simply  a  function  of  the  porosity,  as  defined 
by  Equation  9.  Similar  porosities  in  experiments  4-7  of  about  0.57  give  rise  to  very  different  rates  of 
impregnation.  The  effective  radius  R^,  which  corresponds  to  A  as  calculated  using  Equations  1-3, 
shows  a  much  larger  statKlard  deviation  in  experiments  with  faster  rates. 

Figure  S  shows  the  wicking  curves  for  four  experiments  using  carbon  fiber  tows  having  about  the 
same  porosity.  It  is  possible  to  observe  a  more  gradual  transition  from  linear  to  nonlinear  regimes  of 
impregnation  in  the  faster  experiments.  The  experiments  with  faster  rates  showed  a  second-order  transition 
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Note:  The  bar  denotes  10  urn. 

Rgure4.  Scannine  Electron  Micrograph  of  the  IM7  Carbon  Fibers. 


Table  2.  Analysis  of  the  Data  for  IM7  Carbon  Fibers 


Experiment 

Porosity* 

Liquid” 

Viscosity 

A 

A 

Nonlinear 

Linear 

1 

49.6 

500 

•  0.88 

0.91 

1.55 

0.29 

2 

50.2 

100 

0.89 

0.95 

1.58 

0.33 

3 

51.1 

500 

0.98 

0.99 

1.91 

0.37 

4 

56.8 

100 

1.08 

1.14 

2.32 

0.55 

5 

57.0 

10 

1.46 

1.52 

4.24 

1.25 

6 

57.3 

500 

1.88 

1.88 

7.10 

2.76 

7 

57.5 

10 

0.89 

0.95 

1.58 

0.33 

8 

58.8 

100 

1.16 

1.21 

2.70 

0.49 

9 

59.6 

500 

1.12 

1.07 

2.51 

0.58 

Note:  The  impregnatkm  constants  A  that  are  calculated  using  Equation  19  are  denoted  A  (nonlinear)  and  the  slope  of  the  linear 
portion  of  A  VS.  /  is  A  (linear) 

*  Etpialion  9. 

^  As  dted  by  die  manufacturer. 

*  Weighted  average  of  the  radii  calculated  from  the  non-linear  portion  of  V  VS.  t  (1/n) 

f 

V  * '  ;  _ 

Standard  deviation  of  R^  calculated  from  the  distribution  of  radii  in  a  given  experiment. 
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from  linear  to  nonlinear  behavior  as  compared  with  those  having  slower  rates.  In  the  latter  case,  the 
transition  was  closer  to  being  first-order,  indicating  that  all  portions  of  the  advancing  contact  line  reach 
the  top  of  the  sample  at  roughly  the  same  time.  The  pore  size  distributions  calculated  using  Equation  16 
for  three  of  the  sets  of  data  in  Figure  5  are  given  in  Figure  6.  The  fastest  rate  is  for  experiment  5  where 
the  pore  size  distribution  is  the  widest  and  is  the  largest  These  data  also  allow  us  to  calculate  the 
capillary  number  (Co)  and  thereby  test  our  assertion  that  gravitational  forces  are  negligible  relative  to 
surface  tension  forces.  Letting  Ca  =  2y/prhg,  where  g  is  the  gravitational  acceleration,  we  hnd  that  for 
a  tow  having  a  50%  porosity,  the  average  pore  radius  is  4  •  10'^  mm,  leading  to  Co  =  2  *  10*.  Even  for 
the  largest  pore  diameters  suggested  by  the  distributions  in  Figure  6,  r  =  20  pm  and  Ca  =  2  •  10*. 

When  silicone  fluids  are  used  for  the  impregnation  experiments,  the  initial  linear  relationship  between  V 
and  t  is  valid  for  all  the  carbon  fibers  investigated  in  this  study,  as  seen  from  the  representative  data 
in  Figures  7  and  8.  For  the  epoxy  system,  the  nondimensional  data  corrected  for  non-zero  contact  angle 
also  show  linearity  for  G3O-S00  and  T-IOOOB  fiber  tows  in  the  initial  stages  of  cross-linking  (curing). 
A  slight  curing  effect  is  noticeable  for  both  sets  of  epoxy  data  at  long  times.  The  liiKarity  simply  implies 
that  the  average  cross-set^onal  geometry  in  a  given  tow  is  constant,  although  it  is  nm  possible  to 
determine  the  dqmidetKy  of  the  slope  upon  tiie  geometry  without  a  theory  of  impregnation.  It  is  also  not 
warranted  to  compare  systems  with  different  contact  angles  in  cases  where  normal  impregnation  may  occur 
due  to  the  strong  dependence  of  normal  impregnation  kinetics  on  contact  angle  (Bayramli  and  Powell 
1990). 
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Note:  This  is  assuming  the  Radii  are  discretely  distributed  ova  50  values,  for  three  different  experiments  using  tows  of  similar 
porosities.  The  experiment  niunbers  are  indicated  on  the  curves. 

Figure  6.  Distribution  of  the  Effective  Pore  Radii. 


Note:  Experiments  19  (o).  24(0),  26  (♦)  (0  =  0°),  and  the  epoxy  adhesive  (■),  (6  =  32“).  (Contact  angle  normalization  is 
applied  to  the  epoxy  data.) 

Figure  7.  Linear  Portions  of  the  Impregnation  Curves  for  G30-S00  Carbon  Fiber  Tows.  . 


Note:  Experiments  29  (■),  30  (♦),  (9  =  0“)  and  the  epoxy  adhesive  (O),  (0  =  35°). 

Figure  8.  Linear  Portion  of  the  Impregnation  Carve  of  T-IOOOB  Carbon  Fiber  Tows. 
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The  results  of  the  experiments  with  different  carbon  fibers  are  given  in  Table  3.  The  AS4,  G30-500, 
and  T-1000  B  fibers  are  more  irregular  in  shape  than  the  IM7  fibers.  The  general  trend  of  increasing 
impregnation  rates  with  larger  porosities  is  discernible  from  the  data.  However,  as  in  Table  2,  at  any 
given  porosity,  there  are  many  impregnation  rates.  Experimentally,  we  have  no  control  on  the  porosity 
and  pore  size  distribution  of  the  sample.  Rather,  they  are  measured  during  the  impregnation  experiments. 
Since  the  measurements  of  volume  and,  therefore,  porosity  are  quite  accurate,  the  multiplicity  of  rates  at 
a  given  porosity  is  not  an  experimental  artifact  Rather,  it  represents  a  need  to  characterize  further  the 
interstitial  region  of  the  tow  experimentally  and  to  develop  further  theories  which  will  relate  liquid  uptake 
to  the  pore  geometry. 


Table  3.  The  Slope  of  the  Linear  Portion  of  Impregnation  Curves  for  Different  Carbon  Fiber  Tows 
Having  Various  Porosities 


Run  No. 

CF 

Porosity 

Liquid 

A 

10 

AS4 

45.8 

10 

0.83 

11 

AS4 

53.6 

10 

1.17 

12 

AS4 

56.5 

100 

0.69 

13 

AS4 

57.5 

500 

0.85 

14 

AS4 

60.8 

100 

0.83 

15 

AS4 

64.3 

10 

1.16 

16 

AS4 

64.6 

500 

1.16 

17 

AS4 

65.5 

100 

1.06 

18 

AS4 

69.5 

100 

1.24 

19 

G30-500 

39.6 

100 

0.60 

20 

G30-500 

39.8 

100 

0.47 

21 

G30-500 

41.0 

100 

0.96 

22 

G30-500 

42.9 

10 

0.56 

23 

G30-500 

45.9 

10 

0.52 

24 

G30-500 

49.9 

10 

1.15 

25 

G30-500 

51.4 

500 

0.84 

26 

G30-500 

53.1 

500 

1.02 

27 

G30-500 

55.6 

500 

1.08 

28 

G30-500 

56.2 

500 

1.03 

29 

T-IOOOB 

49.7 

500 

1.07 

30 

T-IOOOB 

56.3 

10 

1.44 

31 

IM7 

38.4 

100 

0.56 
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For  regular  lattices,  the  axial  impregnation  rates  for  different  pore  models  are  shown  in  Figure  9.  The 
experimental  values  of  A  calculated  using  Equation  18  for  porosities  from  0.24  to  0.65  are  also  given  in 
Figure  9.  Equation  8  gives  the  dependency  of  on  porosity,  when  combined  with  Equation  4.  The 
parameters  A,  and  \  correspond  to  the  rate  constants  associated  with  the  pore  radii  r,  and  rj  shown  in 
Figure  1.  They  are  calculated  using  Equations  2  and  3,  with  the  definition  of  porosity.  These 
experimental  vali^s  are  much  higher  than  the  theoretical  ones  based  on  regular  lattice  structures. 
Moreover,  similar  porosities  can  give  rise  to  quite  different  experimental  impregnation  rates.  However, 
for  low  porosities  where  the  packing  density  and,  hence,  the  regularity  of  the  fiber  geometry  is  high,  the 
experimental  rates  are  within  the  limits  predicted  by  the  simple  theoretical  models. 

Cross-sectional  views  of  epoxy-impregnated  fiber  bundles  are  shown  in  Figures  10  and  1 1.  There  are 
regions  of  high  and  low  fiber  density.  There  appears  to  be  no  microscopic  voids  in  the  tow  (a  conclusive 
statement  requires  the  investigation  of  many  cross  sections  of  the  tow;  also,  the  polishing  process  might 
destroy  voids).  The  scanning  electron  micrograph  of  the  tow  cross  section  shows  that  the  fibers  are  indeed 
closely  packed  in  the  high  density  region  (Figure  12).  One  can,  therefore,  postulate  that  the  impregnation 
process  is  governed  by  fast  vertical  wicking  in  the  large  pores  and  a  slow  transverse  penetration  into  the 
smaller  pores  (Bayramli  and  Powell  1990)  which  is  coupled  with  axial  flow  in  these  small  pores.  Such 
a  process  would  give  rise  to  an  apparent  distribution  of  pore  sizes,  as  shown  in  Figure  6. 

5.  CONCLUSION 

These  results  show  that  even  in  the  case  of  impregnation  of  a  simple  porous  medium,  such  as  fibers 
oriented  parallel  to  their  axes,  the  introduction  of  heterogeneity  in  the  fiber-to-fiber  separation  creates  a 
complicated  flow  which  is  very  different  from  that  in  a  perfectly  regular  matrix.  The  parameters  defining 
the  steady  viscous  flow  are  the  porosity,  the  pore  size  distribution,  and  perhaps  the  contact  angle.  Making 
tensiometric  measurements  using  porous  carbon  fiber  tows  allows  the  unambiguous  determination  of 
porosity  and  impregnation  rates  and,  hence,  makes  it  possible  to  elucidate  the  importance  of  these 
parameters  in  the  impregnation  process. 

The  expected  linear  relationship  between  the  volume  of  liquid  impregnated  and  the  square  root  of  the 
time  is  found  for  axially  oriented  fibers  using  liquids  having  contact  angles  of  approximately  zero  and  for 
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Note:  The  numbers  denote  die  experiments  in  Tables  2  and  3;  1M7  (o);  ASA  (>);  G30-S0O  (□};  T-IOOOB  (A).  The  theoretical 
lines  are  for  A,  and  Aj,  csiculated  assuming  cylindrical  capillaries  with  r,  and  r,  (Figure  1),  and  for  A,^  for  the  parallel 
plate  hexagonal  packing  configuration. 


Figure  9.  Comparison  of  the  Theoretical  Dependence 

Lattice  With  the  Experimental]  v  Calculated  Values. 


300  iim 


Note:  The  darkened  boundary  is  the  adhesive  of  the  surrounding  tapie.  The  photograph  is  taken  using  an  optical  reflectance 
microscope.  The  bar  denoted  3(X)  pm. 


Figure  10.  Cross-Sectional  View  of  the  Epoxv-Impreaiated  IM7  Carbon  Fiber  Bundle. 
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Note:  The  phoiogn^  is  taken  using  an  Optical  Reflectance  Microscopy.  The  bar  denotes  200  pm. 


Figure  11.  Cross-Sectional  View  of  the  Epoxv-Impregnated  G30-500  Catfeon  Fiber  Bundle. 


Note:  The  average  fiber  diameter  is  5.0  pm.  The  bar  at  the  bottom  of  the  photograph  denotes  9.7  pm. 

Figure  12.  Scanning  Electron  Micrograph  of  the  Epoxv-lmDregnated  IM7  Carbon  Fiber  Bundle. 


an  uncured  qwxy  system  with  a  finite  contact  angle.  The  observed  rates  of  impregnation  were  sometimes 
over  two  times  higher  than  those  predicted  ftom  the  average  capillary  dimensions  (Figure  9).  Furthermore, 
a  multitude  of  impregnation  rates  is  possible  for  a  given  average  porosity,  differing  by  as  much  as  50% 
in  different  experiments  using  different  tow  samples.  Our  conclusion  is  that  the  average  porosity  alone 
is  not  adequate  to  predict  the  wicking  behavior.  A  knowledge  of  the  pore  size  distribution  is  required  for 
an  accurate  estimation  of  the  impregnation  rates. 

A  large  heterogeneity  in  pore  size  distribution  gives  rise  to  faster  rates  of  impregnation  as  compared 
with  more  homogeneous  distributions.  The  observed  rates  of  impregnation  can  be  qualitatively  explained 
by  a  fast  uptake  of  liquids  in  the  large  pores  coupled  with  lateral  motion  of  the  liquid  into  die  neighboring 
smaller  pores.  It  is  likely  that  a  more  rigorous  analysis  of  the  impregnation  process  can  be  achieved  by 
casting  the  problem  in  terms  of  Taylor-dispersion  relationships  (Brenner  1980)  or  more  sophisticated 
models  of  the  cell  geometry  (Levine  and  Neale  1975).  In  our  case,  the  randomizing  influence  is  the 
convoluted  pore  geometry  in  which  there  is  a  fast  axial  motion  of  the  liquid  in  large  pores  together  with 
an  order  of  magnitude  smaller  motion  in  the  small  pores  of  the  fiber  matrix  as  well  as  lateral  flow.  In 
future  work,  we  will  aim  at  exploring  this  possibility. 

Finally,  the  tensiometric  technique  (measuring  the  amount  of  impregnated  liquid)  and  the  analysis  of 
the  pore  size  distribution  used  in  this  study  can  be  employed  to  characterize  the  structure  of  other  types 
of  porous  material  in  terms  of  wicking  or  wetdng  behavior.  Soil  structures,  solid  porous  catalysts,  or 
woven  fibers  (Kim  1987)  can  be  investigated  with  the  preparation  of  suitable  samples.  Such  smdies  may 
give  clues  as  to  the  dimensional  scaling  of  the  porous  materials. 
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